Thermal properties of a solid through g-deformed algebra 
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We study the thermodynamics of a crystalline solid by applying q-deformed algebras. We based 
part of our study on both Einstein and Debye models, exploring primarily g-deformed thermal and 
electric conductivities as a function of Debye specific heat. The results revealed that q-deformation 
acts as a factor of disorder or impurity, modifying the characteristics of a crystalline structure, as 
for example, in the case of semiconductors. 



I. INTRODUCTION 

Quantum groups and quantum algebras can be con- 
sidered g-deformations of Lie groups and semi-simple Lie 
algebras Physical systems with g-deformation have 
been the subject of intense research due to the emergence 
of quantum group structures caused by certain physical 
problems. This establishes a connection with g-analysis, 
known in mathematics for over a century. From the 
mathematical point of view, g-deformed algebra usually 
requires Hopf algebra. The q-deformation of a physical 
system is not the only way to deform it; rather, there are 
several ways of applying deformation in different physical 
and mathematical contexts 0, Q • 

For example, the widely studied harmonic oscillator 
system has several g-deformed descriptions 0, Q . These 
are obtained from each other by transformation, as shown 
in 0, 0] • Some of the main problems in obtaining the q- 
oscillator are the spectrum, Hamiltonian, position and 
time operator. Moreover, the g-deformation of physical 
systems beyond the oscillator is not well studied. Most 
concepts of classical mechanics and quantum mechanics 
become unclear after being q-deformed Q. 

This study extends our previous analysis [oj to solids. 
A solid consists of a large number of atoms linked by 
cohesive forces of various kinds. On the other hand, in 
a gas, molecules are free to move around the container, 
while in a liquid, they have less freedom, but still travel 
considerable distances. Atomic motion in a solid is very 
slight, causing every atom to move only within a small 
neighborhood, and vibrate around its equilibrium point. 
In a crystalline solid, the equilibrium points of atomic vi- 
brations form a regular spatial structure, such as a cubic 
structure. 

Interaction between atoms allows the propagation of 
elastic waves in solid media, which can be both hori- 
zontal and longitudinal. If atomic oscillations around 
equilibrium positions arc small, which should occur at 
low temperatures, the potential interaction energy can 
be approximated by a quadratic form of atomic displace- 
ments. A crystalline solid, whose atoms interact accord- 
ing to this potential, is called a harmonic solid. In har- 
monic solids, elastic waves are harmonics and the normal 
modes of vibration in crystalline solids [lo| . A large num- 



ber of phenomena involve quantum mechanical motion, 
in particular thermally-activated particles, obeying the 
law. Thermal excitations in the system are responsi- 
ble for phonon excitation (TTj . 

The study conducted by Anderson, Lee and Elliot [l2| - 
[T^ shows that the presence of defects or impurities in a 
crystal changes the electrostatic potential in their neigh- 
borhoods, breaking the translational symmetry of the pe- 
riodic potential. This disturbance can produce electronic 
wave functions located near the impurity, ceasing to be 
propagated throughout the crystal. 

The conductivity of semiconductors can also be dra- 
matically altered by the presence of impurities, i.e., dif- 
ferent from atoms that make up the pure crystal. This 
property enables the manufacture of a variety of elec- 
tronic devices from the same semiconductor material. 
The process of placing impurities of known elements in a 
semiconductor is called doping. Next we see the applica- 
tion of g-deformation in solids acting as a defect or impu- 
rity, specifically in Einstein and Debye solids, and observe 
the consequences; for example, temperature and ther- 
mal conductivity in the Debye model and other features 
of chemical elements We found interesting results, 
leading us to apply the same approach to other thermo- 
dynamic properties of solids, since this initial analysis 
shows that g-deformed elements have thermal properties 
are similar to another element in its pure state. 

The paper is organized as follows: in Sec. |TT]we present 
the g-deformed algebra; in Sec. IIIII we implement q- 
deformation to Einstein and Debye solids in order to ex- 
plore thermal properties; and finally, in Sec. II VI we make 
our final comments. 



II. g-DEFORMED QUANTUM ALGEBRA 

The g-deformed algebra of the quantum oscillator is 
defined by q-deformed Heisenberg algebra in terms of cre- 
ation and annihilation operators and a, respectively, 
and quantum number N by [TB - flTj . 
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where deformation parameter q is real and the observed 
value of q has to satisfy the non-additivity property (see 
3 for a comprehensive study on this property in many 
physical issues) 



[x + y]^ [x] + [y]. 
In addition, the operators obey the relations 
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= [l + iV]. 
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The (/-Fock space spanned by orthornormalized eigen- 
states \n) is constructed according to 



(at) 
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The actions of a, at and N on the states \n) in the g-Fock 
space are known to be 



at|,i) = + + 
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The oscillator 0, [T^ allows us to write the g-deformed 
Hamiltonian [l| as follows: 



n = 



(9) 



We have the basic q-deformed quantum number n defined 
as [11, 



— — = {N\ = a^a = — 



(10) 



At limit g — > 1, the basic q-deformed quantum number 
[x] is reduced to the number x. In our present study 
we shall not attempt to explicitly use Jackson derivative 
{JD) [2 d 1211 to obtain q-deformed thermodynamics re- 
lations [9|, |2^. Instead, we start with the q-deformed 
partition function and use ordinary derivatives to obtain 
q-deformed thermodynamics quantities. 



III. IMPLEMENTATION OF THE 
g-DEFORMATION 

A. g-Deformed Einstein solid 

We consider the solid in contact with a thermal reser- 
voir at temperature T, where we have rij labeling the 
j-th oscillator. Given a microscopic state {uj} = 
{rii, ri,2, . . . , ?T.Ar}, the energy of this state can be written 
as 
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(11) 



Where w^; is the Einstein frequency characteristic. As 
our primary purpose is q-deformation, we proceed to 
find the q-deformed version of (fTT]) . We can obtain q- 
deformed energies from the definition of the Hamiltonian 
(HJ, and the definitions provided earlier, 
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Considering the definition of basic number given in (jlOp , 
and making q = exp(7), we obtain Q 



sinh (1 ) 
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Now replacing the equation ([T^ into equation ((T^), we 
have 



E„ 



sinh 1)1) 
sinh(i) 



The partition function is given by 

S==^cxp(-/3i?K}). 



(14) 
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Since there is no interaction terms, we factored the sum 
and apply the result of equation to obtain the q- 
deformed partition function 
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where 
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As one knows {3 = and Qe is the Einstein tempera- 
ture, defined by 



KB 



(18) 



and we define a q-deformed Einstein function E{aq) as 
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E{aq) = q"' 
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0£(V9+1) _e£;exp(|) + l 



8r^sinh(2) 4rexp(|)-l' 



(19) 



(20) 



We can determine a q-deformed Helmholtz free energy 
per oscillator from the expression 
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By using the result of equation ((2T|) we can determine 
the g-deformed entropy (see Fig. [1]) 
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dT 

^B|ln(<z"')-ln(q"^-l) + ^^|. (22) 



The g-deformed specific heat can now be determined by 
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The specific heat of the Einstein solid as a function of 
E{aq), defined in equation (|f 9p . can be written as follows: 



CVqiT) = 3KBE{aq 



(24) 



The complete behavior is depicted in Fig. [5] One should 
note that when T ^ 9e, with the ratio a = ^ ^ f , with 
9e, for example, around 100 A' for common crystals, one 
recovers the classical result cyq — > 3kb, known as the 
Dulong-Petit law. However, for sufficiently low tempera- 
tures, where T 9e and therefore a 3> 1, specific heat 
is exponential with temperature [ll| as 
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In general, the invariance of specific heat at high temper- 
atures and its decrease at low temperatures show that the 
Einstein model is in agreement with experimental results. 
However, at sufficiently low temperatures, specific heat 
does not experimentally follow the exponential function 
given in equation (j25|) . As for the g-deformed case we 
see a significant change in the curve when g — > 0.1 for 
intermediate temperatures. However, wc will improve 
the model by following the Debye model and applying 
g-deformation in the next section. 

Finally, the g-deformed internal energy per oscillator 
as a function of temperature is given through 
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B. g-Deformed Debye solid 

Corrections of Einstein's model are given by the De- 
bye model, allowing us to integrate from a continuous 
spectrum of frequencies up to the Debye frequency uio^ 



Temperature T(K) 

Figure 1: g-Deformed entropy Sq vs temperature T in the following 
intervals: T = to T = O.IK (top), T = to T = lOOA" (center) 
and T = to T = lOOOE" (bottom). 
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givinff the total number of normal modes of vibration 

mMMM 



g{uj)duj = 



(27) 



where g{u))duj denotes the number of normal modes of 
vibration whose frequency is in the range (a;,a; + duj). 
The function g{ijj), can be given in terms of Rayleigh 
expression as following 



(28) 



where c is the speed of light and A wavelength. The 
expected energy value of the Planck oscillator with fre- 
quency oJs is 



< E, >-- 



1 



(29) 



Using equations ([28| and (|29l). we obtain the energy 
density associated to the frequency range (a;,a; + dw), 
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- 1 
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To obtain the number of photons between oj and w + dw, 
one makes use of the volume of the region on the phase 
space which results in 



g[uj)duj « — 



^''t) (-) 



Vuj^duj 



(31) 



Thus, replacing the equation ([ST]) into equation 
we can write the specific heat for any temperature. We 
now apply g-deformation in the same way as in equation 

(El, 



(32) 



where D(ao)q is the q-deformed Debye function, defined 

by 



"o„ JO 
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Figure 2: g-Deformed specific heat cvq vs temperature T in the 
following intervals: T = to T = O.li^ (top), T = to T = lOOK 
(center) and T = to T = lOOOX (bottom). 
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where uid^ is the g-deformed Debye frequency and is 
the g-deformed Debye temperature. Integrating equation 
(p3| by parts one finds 



D{ao)q 
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exp(aoJ-l al 



a da 



exp(a) — 1 



,(36) 



that can be integrated out to give the full expression 
D{ao)q = 
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exp ( ao„ 
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where 



Li„(z) = 



00 iU 



(37) 



(38) 



is the poly logarithm function. For T S> , then ao, *C 
1, the function D{ao)q can be expressed in a power series 
in ao„ 



D(ao)5 = 1 
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so that for 



(39) 



(40) 



For T <C 0Dg, then ag^ ^ 1, we can write function 
^(ao)^ as 



12 /"^ g^rfg 
ao„ Jo cxp(a) - 1 



+ 0[exp(-aoJ], (41) 
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(42) 



Thus, as in the usual Debye solid, low temperature 
specific heat in a q-deformed Debye solid is proportional 
to rather than proportional to the exponential func- 
tion in ((25|) for the q-deformed Einstein solid. This is 
in agreement with experiments. Finally, we express the 
specific heat for low temperatures as 



= 1944 



T 



J 



molK 



(43) 



For the g-deformed case one can observe the changes that 
occur with Debye temperature, specific heat, thermal and 
electrical conductivies and electric resistivity. By using 
the relationship established for thermal conductivity k 
1251 we obtain 



1 



-CyvL 



(44) 



where v is the average velocity of the particle, Cy is the 
molar heat capacity and / is the space between particles. 
We can deduce a relationship between the thermal and 
electrical a conductivities through the elimination of I (as 
— ne_i ■^yjjQj.f; ^ jg electron mass, n is the number 



FeBnlk CrBulk Feq-Jef - - Ciq-def | 
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Figure 3: Debye temperature Ojj range of Fe and Cr varying as a 
function of g — 0.1 to q 1 (top); deformed specific heat cv range 
of Fe and Cr varying as a function of g — 0.1 to g — 1 (bottom). 
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of electrons per volume unit and e is the electron charge) 
such that 



K 1 CymiP 



(45) 



a 3 ne^ 

In a classical gas the average energy of a particle is 
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Figure 4: Thermal conductivity k range of Fe and Cr varying as 
a function of g = 0.1 to g — )■ 1. 
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kbT, whereas the heat capacity is 



a 3 V e / 
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so 



(46) 



the ratio -% is called the Lorenz number and should be 
a constant, independent of the temperature and of the 
scattering mechanism. This is the famous Wiedemann- 
Franz law, which is often well satisfied experimentally, 
and the Lorenz number correctly given [25|. By using 
the g-deformed relations presented above, we start from 
Eqs. (|44p and P5|) to determine the important relations 
for (jf-deformed thermal and electrical conductivities 



cv 



and 



a, = ^. (47) 



Since electrical resistivity p is the inverse of conductiv- 
ity, g-deformed Pq = Recall that to compute these 
deformed quantities in terms of the specific heat cvq we 
make use of the equation p7p and its suitable limits. 
Data is plotted to provide a better view of our results. 
We see how g-deformation is acting on chemical element 
properties. For illustration purposes, we chose iron (Fe) 
and chromium (Cr) two materials that can be employed 
in many areas of interest. 



Fig. [3] shows how g-dcformation acts on Fe and Cr 
Debye temperatures and specific heat. The plots show 
that Fe reaches rhodium (Rh) values for q k, 0.45, os- 
mium (Os) for q K, 0.25 while Cr approaches silicon (Si) 
for q « 0.42. 

Fig. |4] shows that thermal conductivity for Fe ap- 
proaches the thermal conductivity of rubidium (Rb) for 
q Ri 0.15, whereas Cr equals bulk Fe for q « 0.27, lithium 
(Li) for q w 0.35 and Os for q w 0.48. They also approach 
the germanium (Ge) for q ss 0.17 and q w 0.1, tin (Sn) 
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Figure 5: Electrical conductivity a range of Fe and Cr varying as 
a function of g = 0.1 to q — ^ 1 (top); electrical resistivity p range 
of Fe and Cr varying as a function of 5 = 0.1 to q — > 1 (bottom). 
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Figure 6: g-Deformed thermal conductivity for Cr as 5 — > 1 
and for Cr with impurities as 5 = 0.1 vs temperature T, in 
the interval T = to T = 1500 A' . 



for q K, 0.25 and q «a 0.15, platinum (Pt) and palladium 
(Pd) for q « 0.35 and q ~ 0.18, respectively. 

Fig.[5]shows that electrical conductivity and resistivity 
of Fe equals Ga for q « 0.12, bulk Cr for q « 0.18, Rb 
for q « 0.2, Sn for q w 0.33, Pt and Pd for q « 0.45, 
while Cr equals cesium (Cs) with q « 0.11 and Ca for 
gsa 0.29. 

Fig. |6] shows correct behavior of thermal conductiv- 
ity for a pure and impure material. One should note 
that g-dcformation is clearly playing the role of impurity 
concentration in the material sample. This is because q- 
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deformation acts directly on Debye temperature, which 
means that the Debye frequency is modified. Changing 
the Debye frequency is a clear sign of the material being 
modified by impurities. 



IV. CONCLUSIONS 

Following our previous results in Q, we understand 
g-deformation not only as a mathematical tool, but also 
as an impurity factor in a material, such as disorder or 
reorganization of a crystalline structure. 

In the present study, we first investigate the Einstein 
solid model to obtain thermodynamic quantities such 
as Einstein temperature, Helmholtz free energy, specific 
heat and internal energy with q ^ 1. We then gener- 
alize this study to the Debye solid model to obtain the 
Debye temperature, specific heat, thermal conductivity, 
electrical conductivity and resistivity. We present these 
findings for a number of chemical elements. 

Our main results indicate the possibility of adjusting 
g-deformation to obtain desirable physical effects, such as 
changing the thermal conductivity of a certain element, 
which might become equivalent to a material that is eas- 
ier to handle by inserting an impurity in a sample of the 
original element. We need further studies and evidence 
to substantiate such a complex hypothesis. For example, 
we are seeking to establish a connection between this the- 
ory and experiments through the growth of thin films, a 
matter that will be addressed elsewhere. 
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